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Abstract 

We investigate the spectrum and eigenstates of a Bose-Hubbard chain containing two bosons with fixed boundary 
conditions. In the noninteracting case the eigenstates of the system define a two-dimensional normal-mode space. For 
the interacting case weight functions of the eigenstates are computed by perturbation theory and numerical diagonal- 
ization. We identify paths in the two-dimensional normal-mode space which are rims for the weight functions. The 
decay along and off the rims is algebraic. Intersection of two paths (rims) leads to a local enhancement of the weight 
functions. We analyze nonperturbative effects due to the degeneracies and the formation of two-boson bound states. 
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1. Introduction 

Localization phenomena due to nonlinearity and spa- 
tial discreteness in different physical systems have 
received considerable interest during the past few 
decades. Despite the given translational invariance of 
a lattice, nonlinearity may trap initially localized exci- 
tations. The generic existence and properties of discrete 
breathers - time-periodic and spatially localized solu- 
tions of the underlying classical equations of motion 
- allow us to describe and understand these localiza- 
tion phenomena iH H Si- Discrete breathers were 
observed in many different systems like bond excita- 
tions in molecules, lattice vibrations and spin excita- 
tions in solids, electronic currents in coupled Joseph- 
son junctions, light propagation in interacting optical 
waveguides, cantilever vibrations in micromechanical 
arrays, cold atom dynamics in Bose-Einstein conden- 
sates loaded on optical lattices, among others (for ref- 
erences see m In many cases quantum effects 
are important. Quantum breathers are nearly degener- 
ate many-quanta bound states which, when superposed, 
form a spatially localized excitation with a very long 
time to tunnel from one lattice site to another (for refer- 
ences see |ilj,2i.i4J). 
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The application of the above ideas to normal-mode 
space of a classical nonhnear lattice allowed us to ex- 
plain many facets of the Fermi-Pasta-Ulam (FPU) para- 
dox [5], which consists of the nonequipartition of en- 
ergy among the linear normal modes in a nonlinear 
chain. There, the energy stays trapped in the initially ex- 
cited normal mode with only a few other normal modes 
excited, leading to localization of energy in normal- 
mode space. Recent studies showed that, similar to dis- 
crete breathers, exact time-periodic orbits exist which 
are localized in normal-mo de sp ace. The properties of 
these ^-breathers |S SB 13 [MO, [HQ allow us to 



quantitatively address the observations of the FPU para- 
dox. A hallmark of ^-breathers is the exponential local- 
ization of energy in normal-mode space, with exponents 
depending on control parameters of the system. 

On the quantum side, recently we studied the fate 
of analogous states (quantum q'-breathers) in a one- 
dimensional lattice with two interacting bosons and pe- 
riodic boundary conditions il4ll . By using perturba- 
tion theory, supported by numerical diagonalization, we 
computed weight functions of the eigenstates of the sys- 
tem in the many-body normal-mode space. We did 
find localization of the weight function in normal-mode 
space. However, at variance from the classical case, the 
decay is algebraic instead of exponential. The periodic 
boundary conditions allow us to introduce an irreducible 
Bloch representation. Since states with different wave 
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numbers belong to different Hilbert subspaces, they are 
not coupled by a Hubbard interaction term. Therefore, 
localization along the Bloch wave number is compact. 
This is also happening for the corresponding classical 
nonlinear Schrodinger equation with periodic bound- 
ary conditions lll2ll . when searching for plane- wave-like 
states. 

The classical case however inevitably leads to non- 
compact distributions in normal-mode space, once fixed 
boundary conditions are considered. Indeed, also in the 
quantum case, these conditions violate translational in- 
variance, and lead to nonzero matrix elements between 
states with different Bloch wave numbers, mediated by 
the Hubbard interaction. That is the reason for studying 
the properties of quantum ^-breathers for finite chains 
with fixed boundary conditions. From a technical point 
of view, the irreducible normal-mode space dimension 
is then increased from one to two. 

In Sec. |2]we describe the model and introduce the 
basis to write down the Hamiltonian matrix. We de- 
scribe the quantum states of the lattice containing one 
and two noninteracting bosons. From the latter case we 
use the two-particle states as the basis to write down the 
Hamiltonian matrix in normal-mode space for the in- 
teracting case, after which the energy spectrum is com- 
puted. In Sec. [3]we study localization in normal-mode 
space. We introduce weight functions to describe local- 
ization in that space, and obtain analytical predictions 
using perturbation theory. We present numerical results 
from a diagonalization of the Hamiltonian matrix, and 
compare them with analytical estimates. Then we study 
nonperturbative efi'ects when increasing the interaction 
parameter Finally we present our conclusions in Sec. 

a 



2. Model and spectrum 

We consider a one-dimensional periodic lattice with 
/ sites described by the Bose-Hubbard (BH) model. 
This is a quantum version of the discrete nonlinear 
Schrodinger equation, which has been used to describe 
a great variety of systems 1 15.1 . The BH Hamiltonian is 
H^Ho+yHi Ijil], with 



Ho describes the nearest-neighbor hopping of particles 
(bosons) along the lattice, and Hi the local interac- 
tion between them whose strength is controlled by the 
parameter y. aj and aj are the bosonic creation and 
annihilation operators satisfying the commutation rela- 



tions [3;, fit] 

J 



■'JJ 



[aj,aj,] 



0, and the 



system is subject to fixed boundary conditions. The 
Hamiltonian ([T]l commutes with the number operator 
.^V - whose eigenvalue is n, the total num- 

ber of bosons in the lattice. Here n - 2. It is of interest 
due to its direct relevance to studies and observation of 
two-vibron bound states in molecules and solids 1 18 , 
2fl, 21, 22, 23,24,25,26, 27, 28, 29, 30, 31, 32]. More 
recently, two-boson bound states have been observed in 
Bose-Einstein condensates loaded on an optical lattice 
133]. 

To describe quantum states, we use a number state 
basis |0„) - |ni n2 ■ ■ -rif) where n,- - 0, 1,2 rep- 
resents the number of bosons at the i-th site of the lat- 
tice. lO,,) is an eigenstate of the number operator .^V with 
eigenvalue n - j rij. 

2.1. One-particle states 

For the case of having only one boson in the 
lattice in - 1) a number state has the form 
|0 ■ ■ ■ 1/ ■ ■ ■ 0) = |/>, where I denotes the lattice 
site where the boson is. This number state can be also 
written as 

\i) = a,lo>, (3) 

where the operator a"^ creates a boson at the Z-th site of 
the lattice, and |0) is the vacuum state. 

We write down the Hamiltonian matrix in the ba- 
sis of the above-defined number states. For the single- 
boson case, the interaction term Hi has no contribution 
to the matrix elements. The eigenstates of Hq, for fixed 
boundary conditions, are standing waves: 



/ 



7 + 1 



sin (W) 10 = \k). 



(4) 



where k - qjijif + 1), and q - 1, . . . ,/. The corre- 
sponding eigenenergies are 



Sk — -2cos(A;). 



(5) 



and 



/ 



Hi ^ - Y^a^a^apj. 



(1) 



(2) 



We define bosonic operators aj, fi^ satisfying the 
commutation relations [ak,a^,] - 5k,k', [fl*:,fljt'] - 
[fl^, fl^,] = 0, such that the state (|4]i may be written sim- 
ilar to ©: 



.ka, 



(6) 



z=i 



where the operator creates a boson in the single- 
particle state with quantum number (wave number or 
momentum) k. The bosonic operators at, are related 
to the operators S;, in direct space through the trans- 
formation matrix 



sin(fc/). 



(7) 



2.2. Two-particle states 

For the two-boson case (n - 2), we define the number 
state basis in a similar way as in the single-boson case: 



2 — 5/| /, 

^ a;,a/,|0), 



(8) 



where l2> h because of the indistinguishability of par- 
ticles, a^^ and a^^ respectively create one boson at the 
lattice sites h and The number of basis states is 
d = f(f + l)/2. The interaction term Hi in ([TJ con- 
tributes to the matrix elements of the Hamiltonian in the 
above-defined basis. 

In the noninteracting case (y - 0) the eigenstates of H 
in terms of bosonic operators in the normal-mode space 
read [see Eq. ©I: 



1^1,^2) 



2-6, 



- «i,«fclO) , qi>q\ 



(9) 



a^^ and a^^ respectively create one boson in the single- 
particle states ki and k2 of the form Using Eqs. Q 
and (|7|, the relation between the basis states in normal- 
mode space (|9]l and the basis states in direct space © 
reads: 



\kuk2) 



V2 - ^^1,^2 

V2 

/ / 

h = lh>li 

f 1 

+ V2^5,,,,5yJ/,/) . (10) 

1=1 

In the interacting case (y > 0), we represent the eigen- 
states of the Hamiltonian ([TJ in the normal-mode ba- 
sis ( fTOb of the noninteracting case. This leads to a 
d X d matrix [d - f(f + l)/2] whose elements H{i,j) 
1,. . .,d) are 

H(i,j) = {k\,k'^\H\kuk2) = {q[,q'2\H\quq2). (11) 

The integer j that labels the column of the matrix ele- 
ment (fTTT i is related to the mode numbers qi and q2 by 



y?i,?2 -!)(/+ !)■ 



(^1 - l)(qi+2) 



+ q2. (12) 



The same relation holds for the integer labeling the 
row of the matrix element ( fTTT i. 
The matrix elements (fTTT i are 



where 
and 



H{i,j) = HQ(iJ) + jHiii,j), 



Ho(i,j) = (sk, + Sk^)6ij, 



(13) 



(14) 



Hid, j) - fqi,q2,<]\,l'2 ^ ^ IM^Ukl^ l,k\S l,k',- (15) 



1=1 



Sk is the single-particle energy given by Eq. and the 

coefiicients/^,,^,v,f/ are 



fll,l2 



9l>92,9,,92 



^(2-5,„,,)(2-V„,;) 
(/+1)2 ■ 



(16) 



In Fig. [T]we show the energy spectrum of the Hamil- 
tonian matrix ( fTSl l obtained by numerical diagonaliza- 
tion for dififerent values of the interaction parameter y. 
In all calculations by numerical diagonalization we used 
/ = 40, which leads to a matrix dimension d = 820. The 
eigenstates are ordered with respect to their eigenvalues 
Ey (v - I, . . .,d). At y = 0, the spectrum consists of 
the two-boson continuum, whose eigenstates 1^1,^2) are 
given by ( fTOb . The eigenenergies are the sum of the two 
single-particle energies: 



^k,,k2 



- -2[cos(ki) + cos(k2)]. 



(17) 
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Figure 1: Energy spectrum of the two-boson BH chain with fixed 
boundary conditions for different values of the interaction strength y. 
The eigenvalues are plotted as a function of the eigenvalue label (see 
text). Here / = 40. 



When y > 0, eigenvalues in the lower part of the spec- 
trum are pushed down, and beyond y » 2 a band of / 
states splits off from the two-boson continuum. These 
are the two-boson bound states, with a high probability 
of finding the two bosons on the same lattice site, while 
the probability of them being separated by a distance r 
decreases exponentially with increasing r 

The critical value yi, = 2 for which the band of two- 
boson bound states splits off from the continuum may 
be explained as follows. In the limit / oo the un- 
normaUzed bound state with highest energy E = -2y is 
given by ifllflTll: 



l*P> = Xi(-l)UO. 



(18) 



For = 2 the energy of that state leaves the two-boson 
continuum of energies E e [-4, 4]. 



Thus for ki + k\ and k2 + ka the weight function 
C(ku k2; kuh) ^\{ku kim^ljl^ is 



C(ki,k2;ki,k2) - y' 



\{kuk2mh,~k2}\^ 



(21) 



where ii^^ and are eigenenergies of the unper- 



turbed system given by ( fTTl i. For convenience we use 
new variables in normal-mode space 



kj, = k2 + ki, 



(22) 



which are the total (Bloch) and relative wave numbers 
respectively. They have values < k+ < In and < 
k- < n. Since we are interested in the behavior of the 
weight function around the core at (^i , ^2), we define the 
coordinates relative to that point: 



(23) 



3. Localization in normal-mode space 

We recall that the normal-mode space is spanned by 
both momenta ki and k2. The conditions < ^i_2 < 
and k\ < k2 reduce the normal-mode space to a triangle 
that we call the irreducible triangle, as sketched in Fig. 
|2l For finite / and y the eigenstates |*P) will spread in 
the basis of the y = eigenstates A:2)). We measure 
such a spreading by computing the weight function in 
normal-mode space C{ki,k2) - \{ki,k2^)'^ . 

3.1. Analysis by perturbation theory 

We use perturbation theory to calculate the weight 
functions, where y is the perturbation. We fix the mo- 
mentum ki and k2, and choose an eigenstate 1^1,^2) of 
the unperturbed case y - Q. The wave numbers ki and 
k2 define a seed point P - (ki, ^2) in the irreducible tri- 
angle (see Fig. |2]i. Upon increase of y, the chosen eigen- 
state transforms into a new eigenstate l^l'jjj,), which will 
have overlap with several eigenstates of the y = case. 
We expand the eigenfunction of the perturbed system to 
first order in y: 



Thus, (I2TI) becomes 



l^iJ2) + yl*Pf. ), 



(19) 



where 



k'>k', 



{k\,k'mkuk2) 



£9_ 

k,k2 



k'.k' 



\k\X). (20) 



C(k\,k2\kx,k2) - Y 



2 •'91,92,91,92 



[16(£° -E°)]^ 

^ k,k2 hk2'^ 



X<,._;^.,L' (24) 

where /^,,, 2,91, 92 is given by Eq. (O. 

The coefficient 7?^.^ ^ i_ consists of a sum of eight 
terms of the form 

sin[(2/+l)fl 

giO = ^ , (25) 

sin (I) 

with pairwise opposite signs (see appendix |A|. For 
each term, the argument f is a certain combination of 
the wave numbers k+,k- and k+,k^ (see appendix lAl 
for details). Unless the argument of any of the eight 
terms g(^) vanishes, all of them cancel each other and 
^k^ k_-k+ L - 0- Thus the condition ^ = for each 
term in 7?^ ^ .j j , together with the relations (l22l l and 
(l23T l. defines lines k2 - ^2(^1) in the normal-mode space 
where the weight function C{ki,k2\k\,k2) is nonzero. 
These lines are schematically shown in Fig. |2](the an- 
alytical derivation of these lines is given in appendix 
IaI . Note that these lines are specularly reflected at the 
boundaries = and k2-nof the irreducible triangle. 

To study the localization in normal-mode space away 
from the core using the formula (|24] |. we consider the 
two cases A_ = 0, A+ > and vice versa, i.e. the mu- 
tually perpendicular directions A+ and A_ (Fig. |2]i. For 




Figure 2: Sketch of the different lines in the two-dimensional normal- 
mode space along which the weight function 1211 is nonzero. The seed 
point P = (k\,k2) corresponding to the unperturbed eigenstate \k\ , ^2) 
is represented by the black spot. Its conjugate point P = {7T-k2,n-k\) 
is represented by the grey spot. The axes defining the coordinates A-,- 
and A- are indicated by the arrows emerging from P. 



each case we obtain, with |A±| «c n. 
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AJI [cos(fei) -i- cosfe)] 



-1- sin(A:i) + sin(A;2) 



(26) 



The effective interaction strength is y/if + 1). In the 
hmit J or / ^ oo we have compactification of 
the eigenstates. The formula ( |26] | shows locaHzation in 
normal-mode space. Depending on the seed (^i , ^2) we 
find algebraic decay within the irreducible triangle, C ~ 
A-", with a = 2,4. If sin^i ± sinA:2 0, a = 2. If 
sin^i ± sin ^2 -0, a = 4. E.g. for^i = ^2 



C 



J 



1 



/+ 1/ cos2(fei)A;' 



(27) 



Note that along the A+ direction in the irreducible tri- 
angle, Rk^^tSk+,k- ~ 2(/ + 1) at all points but P - (k\ - 
7T -k2,k2 - 71 - ki). This is the conjugate point of the 
seed P (Fig. |2|, where two lines intersect. At this point 
k_-k+ k_ - 4(/ + !)■ Thus we expect a local maximum 
of the weight function at the conjugate point. The states 
1^1, ^2) and 1^1, ^2) have energies E' 







kuh 



3.2. Numerical results 

In Fig. |3] we show the weight function in the two- 
dimensional normal-mode space obtained by numerical 



Figure 3: 3-D plot of the logarithm of the weight function in the 
normal-mode space for the eigenstate v = 145, obtained by (a) nu- 
merical diagonalization, and (b) perturbation theory using the formula 
f24l . Here / = 40 and y = 0.1. 



diagonalization and the formula (|24] | respectively, with 
characteristic localization profiles. We find agreement 
of the numerical data with the results from perturbation 
theory. The largest value is at the point P - ( jqTi, ^n) ~ 
(0.27r, 0.47r), and it decays mainly along the lines de- 
scribed in the previous section (Fig. |2]). Note also the 
presence of the local maximum at the conjugate point 
P ~ (Q.6n, O.Stt) in both cases. 

In Figs. |4]and|5]we plot the weight function of the 
eigenstate shown in Fig. |3] along the directions A+ and 
A_ respectively for different values of the interaction pa- 
rameter y. The state becomes less localized with in- 
creasing y, as expected from the above analysis. The 
decay of the weight function is well described by per- 
turbation theory (dashed lines). The peak of the weight 
function at the conjugate point is clearly seen in Fig. HI 

In Fig. |6] we plot the weight function of different 
states along the A+ direction. It decays as a power law 
that ranges from A""^ for states near the lower corner of 
the irreducible triangle (see Fig. [8]) to A"^ for states ful- 
filling k2 ~ n - ki. In Fig. |7]we plot the decay of the 
weight function along the A_ direction, where we see 
the power-law decay that ranges from A"^ for states ful- 
filling ki X k2 (see Fig. [8|, to A~^ for states fulfilling 




Figure 4: Weight function for different values of the interaction 
strength y of the eigenstate v = 145 along the A+ direction. The 
dashed lines are results from formula 1241 . Here / = 40. 
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Figure 7: Weight function of different eigenstates (labeled by the in- 
dex v) along the A_ direction. Here y = 0. 1 and / = 40. 
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Figure 5: Weight function for different values of the interaction 
strength y of the eigenstate v = 145 along the A_ direction. The 
dashed lines are results from formula 1241 . Here / = 40. 
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Figure 6: Weight function of different eigenstates (labeled by the in- 
dex v) along the A+ direction. Here y = 0. 1 and / = 40. 



k2 ~ 7T-k[. The results from numerical diagonalization 
agree very well with those from the perturbation theory 
analysis. 



Figure 8: Location P = (^1,^2) of the eigenstates, shown in Figs. |6] 
and0 in the irreducible triangle. 



3.3. Nonperturbative effects 

The results in the previous section were obtained 
for small values of the interaction parameter y up to 
■y - 0.1, for which perturbation theory gives a good 
description of the results obtained by numerical diag- 
onalization. However, when increasing y several non- 
perturbative effects occur These are: 

Split off of the two-boson bound state band: This ef- 



feet was discussed in Sec. 12.21 (Figs. I 111. When y > 2 
the two-boson bound state band splits off from the two- 
boson continuum, and the corresponding eigenstates are 
correlated in direct space, i.e. with large probability the 
two bosons are occupying identical lattice sites. Thus, 
in normal-mode space these eigenstates become delo- 
calized as shown in Fig. |9] 

Degenerate levels in the noninte meting case: The 
analysis using perturbation theory is valid as long as the 
eigenstate which is continued from the noninteracting 
case is not degenerate. Because of the finiteness of the 
lattice the momenta ^1 and ^2 are restricted to discrete 
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Figure 9: 3-D plot of the logarithm of the weight function in normal- 
mode space for the eigenstate v = 25, that belongs to the two-boson 
bound state band. Results were obtained by numeiical diagonaliza- 
tion. (a) 7 = 1.5. The two-boson bound state band did not split off 
and the eigenstate is localized in normal-mode space, (b) y = 2. At 
this interaction value the two-boson bound state band splits off and the 
eigenstate becomes delocalized in normal-mode space. Here / = 40. 



values and define a grid in the two-dimensional normal- 
mode space. A grid point (^i, ^2) defines a line of con- 
stant energy in normal-mode space through Eq. (fTTI l. 
with = £? - (Fig. [TOla). The nondegeneracy 

condition implies that this line should not pass through 
any other grid point. It is easy to see from Eq. (flTl l that 
all states |^i,7r - ^i) are degenerate, with - 0. 

Their corresponding grid points in the irreducible trian- 
gle lie on the diagonal k2 - n - k\ (thick line in Fig. 
[TOl-a). In Fig. [TOl -b we show the weight function of an 
eigenstate that is located on that diagonal in the non- 
interacting case. As expected, even for small values of 
7, the state completely delocalizes along the degeneracy 
diagonal. 

Avoided crossings: Upon increase of the interac- 



tion parameter y, the energies of continued eigenstates 
change, and will resonate with eigenvalues of other 
states. 

The first possible avoided level crossing defines a 
critical value of the interaction parameter y up to which 
first-order perturbation theory is applicable. To estimate 



Figure 10: (a) Lines of constant energy given by U7t . The red thick 
line is the i;^^ ^.^ = 0-line where all unperturbed states \ki,k2 = n-k{) 
are degenerate, (b) 3-D plot of the logarithm of the weight function in 
normal-mode space for the eigenstate v = 405, obtained by numerical 
diagonalization with 7 = 0.1 and / = 40. In the noninteracting case 
this state corresponds to\k\,k2 = n - k\). 



this value, yc, we assume that before the first avoided 
crossing is encountered, the eigenenergies depend lin- 
early on y. This dependence may be estimated using 
first-order perturbation theory in y. The result is, for 
large /, 

where 

( 2 ifli = 0, 

b{ku~k2)^\ -2 if^2=^i>0, (29) 

{ -1 if > 0, ^2 > ^1- 

Let us consider two levels £1 and E2 that interact in 
the first avoided level crossing. At 7 = they are sepa- 
rated by 6E. For nonzero y the energies linearly change 
in y: 



\h + 7 ^' 



17 I7O , cc- 1^2(^1,^2)1 

^2 = ^i„fc + ^^ 7 — y- 



(30) 



(31) 



By equating £1 and Eoatj - jc we obtain 

f6E 



7cik\,k2) = 



\bl(kuk2)\ + \b2(kuk2)\ 



(32) 



The first avoided crossing of the level £1 will happen 
with its nearest neighbor in the spectrum of the y = 
case, which is separated by 6E(k\,k2)- Using Eq. ([TtT i 
with k\,k2 ^ 7r/2 (see Fig. [TOl-a). the separation is esti- 
mated (see appendix iBt: 5E k 4-"^n/(f +1)^. There- 
fore 7^-!//- 

The coefficient bi(ki,k2) depends on the state ^2) 
under consideration through ( |29] l. The coefficient 
^2(^1,^2) must have opposite sign as compared to 
t>i(ki,k2) for the avoided crossing to take place. For 
the states v = 145 and v - 41 located at (^1,^2) ~ 
(0.27r, 0.47r) and (^1,^2) ~ (0.12;r, 0.227r) respectively 
(see Fig. [8|, bxikuh) = -1 and ^72(^1, ^2) = 2. This 
leads to a critical value of the interaction parameter 
Jc ~ 0.28, which is in reasonable agreement with the 
numerical results: jc ~ 0.2 for the state v — 145, and 
Jc ~ 0.3 for the state v = 41. 

4. Conclusions 

In this work we studied the properties of quan- 
tum q-breathers in a one-dimensional lattice contain- 
ing two bosons modeled by the BH Hamiltonian with 
fixed boundary conditions. Because of the lack of 
translational invariance, the normal-mode space is two- 
dimensional and reduces to a triangle when working in 
the irreducible representation of the product basis states 
(the irreducible triangle). To explore localization phe- 
nomena in this system we computed appropriate weight 
functions of the eigenstates in the normal-mode space 
using both perturbation theory and numerical diagonal- 
ization. We find that the weight function is sizable only 
along the mutually perpendicular directions defined by 
the total and relative momentum, thus it defines lines in 
the irreducible triangle that show specular reflections at 
the boundaries of the irreducible triangle. We observe 
localization of the weight function along these lines. 
The localization is stronger when the size of the sys- 
tem increases or the interaction parameter is weaker, the 
former because the effective interaction drops in the di- 
lute limit of large chains. We found algebraic localiza- 
tion. The power of the decay is different for each eigen- 
state depending on which seed wave numbers have in 
the noninteracting case, ranging from two to four 

An interesting effect is the local maximum of the 
weight function at the symmetry-related (conjugate) 
point of the eigenstate core in normal-mode space, due 



to a crossing between different paths described by the 
lines along which the weight function is nonzero within 
perturbation theory. 

In addition to the existence of degeneracies between 
eigenstates in the noninteracting case, we analyzed 
other nonperturbative effects as the interaction param- 
eter increases, which limit the applicability of perturba- 
tion theory to describe the system: The splitting off of 
the two-boson bound states from the two-boson contin- 
uum, and the occurrence of avoided level crossings. The 
first effect manifests as a delocalization of the weight 
function of the bound states due to the two-boson cor- 
relation in direct space. The second effect manifests as 
a sudden change of the location of an eigenstate in the 
normal-mode space due to resonant interaction with an- 
other eigenstate. Both effects define critical values of 
the interaction parameter below which one may analyze 
the system by perturbation theory. The occurrence of an 
avoided level crossing gives the smallest critical value. 

Although we considered a system with fixed bound- 
ary conditions, we still obtain algebraic decay as in 
the case with periodic boundary conditions lll4ll . The 
question how to restore exponential localization of clas- 
sical q-breathers from algebraic decay of quantum q- 
breathers in the limit of large numbers of particles is still 
open. WTien going to that limit, one may use a Hartree 
approximation and describe the system with a product 
state wavefunction, or use a coherent state representa- 
tion. Both ways lead to the nonlinear Schrodinger equa- 
tion where classical q-breathers are known to exist lll2[l . 
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A. Lines of nonzero weight function 

For fixed k\,k2, the coefficient R(k\ ,k2',k\, ^2) in Eq. 
(l24l l is given by 

RihMMM) = g{k^+kS) + g{A_) 

- g{k+ + L) - g(k+ - L) 

+ g(A^)+g(k^+hX (33) 

where 



8(0 



sin[(2/+l)|] 
sin(f) 



(34) 



The lines k2 = ^2(^1) in the normal-mode space (irre- 
ducible triangle) along which R{ki,k2\k\,k2) are 
obtained from the condition that the argument of any 
term in (l33T l is zero, such that 

g(0 = 2/+l. (35) 

Let us analyze each of the arguments in Eq. (l33T l: 

• ^_ + ^_ = 0: This implies that 

k2-ki=-(k2-~ki). (36) 

Since k2 > ki the above condition is possible only 
for points (^1 , ^2), (^1 , ^2) on the diagonal ^2 - ki. 



A_ = 0: This condition leads to 

k2 ^(k2-ki) + ku 



(37) 



which is the equation of the line along the A+ di- 
rection that cuts the ^2 axis at k2(Q) -k2-k\. 



k- +k+ - 0: This implies that 

k2-kx - -(k2 + ki), 



(38) 



which is possible only if ^1 = ^2 = and (ki , ^2) is 
on the diagonal k2-k\. 



k- -k+ - 0: This leads to the equation 

h ^{k2 + h) + ki. 



(39) 



which describes a line parallel to the A+ direction 
that cuts the ^2 axis at A;2(0) -k2+k\. 

k+ +k^ - 0: This implies that 

k2^-{k2 + h)-ki, (40) 
which is valid only if k\ - k2-k\ -k2- 0. 
k+-k-- 0: This leads to the equation 

k2^(k2-ki)-ku (41) 

which is the equation of a line parallel to the A_ 
direction that cuts the ^2 axis at k2{Q) -k2-k\. 



A+ = 0: This leads to the equation 

k2 ^{k2 + ~ki)-ku 



(42) 



which describes the line along to the A_ direction 
that cuts the ^2 axis at ^2(0) -k2+ki. 
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Figure 11: Sketch of the discrete normal-mode space (irreducible 
triangle) with a line of constant energy in the circular approxima- 
tion passing through the grid point P = (k[,k2) and the grid point 
(ki + Ak, k2 + ^k). The strip of width d = ijlhk and area A , con- 
tains Ng = A JAk^ grid points through which Hnes of constant en- 
ergy pass. The typical line separation within the strip is Sk x: d/Ng. 
Ak = !T/{f + 1) is the grid spacing. 



B. Energy separation between nearest-neighbor lev- 
els in the noninteracting case 

The finite size of the lattice leads to discrete values of 
the momenta ki and k2, and thus to a grid in the normal- 
mode space (Fig. fTTI) . Let us consider a line of constant 
energy which passes through the seed point P - {k\, ^2) 
given by £? , = r with Eq. ( fTTI ). For small values 
of k\ and k2, the energy in Eq. (fTTb may be approxi- 
mated to 

(43) 



h,k2 



-A + k\ + kl, 



which is the equation for a circle (circular approxima- 
tion). So the equation for the line of constant energy 
passing through the point P is 



k2{ki;k\,k2) 



■kl. 



(44) 



Through another grid point at (ki + Ak, k2 + Ak), sep- 
arated from P by a distance d - y/lAk ^ Ak (Ak is the 
grid spacing), another line of constant energy with the 
form ( l44l i passes (Fig. fTTT i. defining a strip of area Aj 
in the irreducible triangle. The strip contains Ng grid 
points through which lines of constant energy pass. The 
average line separation within the strip is 6k x d/Ng. 

The number of grid points in the strip is 



N, = 



Ak^' 



(45) 



The area of the strip is 



;[(ki + AkY + ik2 + Akf -k]- kj] 

I 

'-[Ak^ + (ki +k2)Ak]. 



Therefore 



and hence 



n 

4Ak' 



Ng - -Y^{Ak + ki+ k2). 



4Ak^ 



Sk(kuk2) = 

7T{ki + k2 + Ak) 

The corresponding energy separation is 



(46) 
(47) 

(48) 



^^(^iJ2) = £°,,,„,^,,,^-4„,^, (49) 

with 6ki - 6k2 - 6k/ V2. Substituting (|43]l and (l48T l 
into ( |49l ) one obtains 



6E(kuk2) 



16Ar 



TiHki +k2 + AkY 

4 V2(^i + k2)Ak^ 
n(k\ +k2 + Ak) 



(50) 



For / large, and Ak - njif +\) <^ fei,2, the first term in 
I can be neglected. Thus we are left with 



5E{ki,k2) 



(/+1)2 



(51) 
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